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The Higgs boson may dominantly decay to 4 light jets through a light pseudo-scalar intermediary: 
/i — >■ 2r; — >■ 4j, making reconstruction at the LHC particularly challenging. We explore the phe- 
nomenology of such "Buried Higgs" scenarios in which the primary discovery channel of the Higgs 
is in cascade decays of superpartners. QCD backgrounds that would otherwise overwhelm the Higgs 
decay are suppressed by the requirement of high pr jets and large missing transverse momentum 
that are the typical signatures of TeV scale supersymmetry. Utilizing jet substructure techniques, 
we find that for buried Higgses in the 100 — 120 GeV range, a 5a discovery can be expected with 
roughly 10 — 25 fb~^ of data at Ecm = 14 TeV. For lighter Higgs bosons, the signal is contaminated 
by hadronically decaying W bosons, and discovery remains an unsolved challenge. 



I. INTRODUCTION 

Understanding the mechanism of electroweak symme- 
try breaking (EWSB) is the primary goal of the Large 
Hadron CoUider (LHC) program. Within the Standard 
Model (SM) of particle physics, EWSB occurs via the 
Higgs mechanism which predicts a weak scale Higgs bo- 
son h which unitarizes longitudinal W boson scatter- 
ing and accommodates electroweak (EW) precision data. 
The current limits on the SM Higgs boson mass rrih con- 
sist of the lower bound nifi > 114 GeV from LEPH [1], 
and the exclusion window 158 < rrih < 175 GeV from 
Tevatron [5]. However, these limits are model dependent. 
For instance, in several extensions of the SM with an ex- 
tended Higgs sector, new decay channels for the Higgs bo- 
son are available and the exclusion bounds are weakened 
[3H5] . Therefore, new strategies to discover non-standard 
Higgs bosons at the LHC and the Tevatron need to be de- 
veloped to be sure all regions of EWSB parameter space 
are covered. 

In this work wc consider the possibility that a light 
Higgs boson is produced in supersymmetric cascades and 
subsequently decays to four light jets (i.e without 6's), 
/i — > 2ry — )■ 4j, where ?7 is a new light pseudoscalar that 
is neutral under the SM gauge group. This non-standard 
Higgs boson decay is present in many extensions of the 
SM that stabihze the EW scale [3 [7HT^. The typical 
fine-tuning in these models is reduced because the exper- 
imental constraints on the Higgs sector are much weaker. 
Indeed, for 77 heavier than 12 GeV only the OPAL model 
independent bound [13], ruh > 82 GeV, holds. If < 12 
GeV, the lower bound is slightly stronger [14], m/i > 86 
GeV, while recently recasted [H] ALEPH data [16] set 
rUrj > 4 GeV for rrih — 100 GeV. Present hadron collider 
data do not constrain such a Higgs boson because of the 
huge QCD background. 

Recent studies [TT] [TH] have discussed search strategies 



of a Higgs boson decaying to four light jets at the LHC by 
means of jets substructure techniques [T^]. ^ If the Higgs 
boson has a large boost factor, its decay products will be 
merged into a single fat jet with characteristic substruc- 
ture. Since QCD jets from background processes are rela- 
tively uncorrelated, substructure requirements can signif- 
icantly enhance the signal to background. Our analysis 
extends previous works in three main directions: 

• We consider Higgs boson production in association 
with supersymmetric particles which can lead to a 
much earlier discovery (similar to the h ^ bb anal- 
ysis in SUSY [211 122] ) compared to the search using 
SM-like production. In particular, we consider neu- 
tralino and chargino two body decays, Xi ~^ x'jh 
and -^'xfh. 

• We do not restrict our analysis to very light 77 below 
the bb threshold. We instead extend the -q mass 
range, m,, < 30 GeV. 

• We discuss how to reduce the contamination from 
hadronic W and Z for a low-mass Higgs boson. 
These are the irreducible backgrounds for super- 
symmetric Higgs productions and can be larger 
than the signal in the resonance regions. This re- 
quires special treatment because they can easily 
produce a similar substructure pattern as the Higgs 
boson. 

We find that for rrih < 95 GeV and for a light pseudo- 
scalar = 10 GeV, the Higgs boson remains hidden in 
the W background after all cuts are applied. For heavier 
Higgs bosons, where the VF's are well separated from the 



For an alternative search based on tagged outgoing forward pro- 
tons see I20| . 
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Higgs boson peak, we obtain 5a signal significance at 
the LHC with = 14 TeV and 10 - 25 fb"^ which is 
smaller by a factor 5 — 10 than the luminosity needed for 
discovery in the SM production channel \TP, 15] . 

This paper is organized as follows. In Sec. |ll] we dis- 
cuss some illustrative models in which the Higgs origi- 
nates in a supersymmetry cascade and subsequently de- 



cays to four light jets. In Sec. |IIH we discuss rates of 
Higgs boson production associated with supersymmetric 
particles. These supersymmetric events are isolated by 
means of cuts discussed in Sec. IIVI In Sec. |V] we discuss 
our analysis, including a refinement of jet substructure 
algorithms used in our search, and we present the results 
in Section [¥!} Finally, we conclude in Sec. |VII| 



II. HIGGS DECAY TO FOUR JETS 

Cascade decays of the Higgs boson to SM states 
through a pair of light pseudoscalar 77's are well moti- 
vated in several extensions of the SM. They include su- 
persymmetric [SirfMllj as well as non-supersymmetic |12j 
realizations. After spontaneous EWSB, a trilinear cou- 
pling between light SM singlets rj and the Higgs boson h 
is a generic feature: 



(1) 



where / is the typical scale that controls the interaction 
strength. For instance, in models where h and 77 emerge 
as pseudo-Goldstone bosons and <C mh, the deriva- 
tive interaction h{d^rj)'^ahriril f produces the Lagrangian 
([1]) after integration by parts. In these scenarios / is anal- 
ogous to the pion decay constant that sets the couplings 
in the chiral Lagrangian. 

If / is not too large compared to the EW scale and 
afiriri is not too small, the decay width of the Higgs boson 
into two 77's 
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can easily dominate over the SM channel h — )■ bb. The 
pseudoscalar rj is generically unstable because it cou- 
ples to SM fermions with effective Yukawa interactions 
iytprj^p'-f^'ip generated after EWSB, « m^/^/2f. Then, 
barring accidental cancellations, the largest Yukawa cou- 
pling is to the third generation quarks, while all other 
SM fermion couplings are highly suppressed. For re- 
cent studies where 77 — )■ &6 see Ref. However, 
when rrirf < 2mi, ^ 10 GeV the decay of 77 to two 
gluons via loops of third generation quarks will be the 
dominant decay mode, corresponding to a four unfla- 
vored jet final state, /i — ?> 277 — ?> Ag. For instance, 
this scenario is naturally realized in the supersymmet- 
ric "Buried Higgs" model [7 where both h and 77 are 
pseudo-Goldstone bosons arising from a global SU{3)h 



symmetry broken down to SU{2)h at the scale / « 500 
GeV. In this model, the coupling between r? and h de- 
pends only on a mixing angle v/ f which measures the 
alignment between the gauged SU{2)w and the residual 
global SU{2)h. aunr, « v/{V2J){l - / P)-^'\ The 
branching ratios for ft, — > 277 and 77 — > 2^ arc 80 — 90% 
and 100% respectively. 

An appealing modification of the Buried Higgs scenario 
is the "Charming Higgs" model 0. The Higgs sector is 
the same as in the original Buried Higgs model but the 
embedding of the matter content into the S'[/(3)h global 
symmetry multiplets is different. In particular, the bot- 
tom Yukawa arises only from non-renormalizable opera- 
tors suppressed by the physical scale A « 10 TeV where 
new heavy states are integrated out. Then, the resulting 
bottom Yukawa coupling to 77 is greatly suppressed be- 
cause it has to vanish both for large / and for large A, 
therefore y?, ~ mh/ ^/2f xmH <C 1 . Thus, it turns out 
that the dominant decay channel is the tree-level 7/ — >■ cc 
even when 77 is above bb threshold production, > 2m;,. 
The next relevant decay mode, 77 — 2g, is generated at 
1-loop. Very much like the original Buried Higgs model, 
the charming version buries the Higgs boson beneath the 
QCD background at the LHC. 

Another class of models where the Higgs may naturally 
cascade decay to four jets is provided by the non-minimal 
composite Higgs models iX2\. 

In the following, we will use the Buried/Charming 
Higgs models as illustrative examples of supersymmet- 
ric extensions of the SM where the Higgs boson cascade 
decays to four light jets where € [5 — 30] GeV and 
777,;i e [90 — 120] GeV. For simplicity and clarity of pre- 
sentation, in our simulations we consider Higgs boson 
production rates as they appear in the Minimal Super- 
symmetric Standard Model (MSSM). ^ In particular the 
Higgs production from SUSY events is matched to the 
production in the MSSM with the same input parame- 
ters and Higgs mass. In our numerical study we take the 
branching ratios for h ^ 21] and 7/ — )• 2j to both be 100%. 



III. HIGGS PRODUCTION 

In general, there are several possible production chan- 
nels for the Higgs boson in SUSY in addition to the 
SM channels. The important ones are those with large 
cross sections such as the pair production of gluinos and 
squarks. In the subsequent cascade decay of these par- 
ticles, a Higgs boson can be produced in many different 
ways, the most important being through supersymmet- 
ric gauge interactions from the gaugino-Higgsino-Higgs 
coupling. Depending on the mass spectra of charginos 



^ Up to / corrections, these rates are the same as those in the 
buried and charming Higgs models. 
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and neutralinos, there can be two scenarios: the "httle 
cascade" and the "big cascade" . 

In the first Higgs is produced in the decay 

Xi ~^ Xih- For this channel to dominate, the /i-term 
needs to be larger than the gaugino masses to ensure that 
Xi and Xi 2 ^-^'S gaugino-like. If, in addition, the sleptons 
are heavier than the decay to Higgs is always domi- 
nant since the other available mode, X2 ~^ Xi^ depends 
on the Higgsino-Higgsino-Z coupling which is doubly sup- 
pressed by small Higgsino mixing. The X2 can be pro- 
duced from the decays of g or g with sizable branching 
fractions, producing a sizable fraction of SUSY events 
which contain Higgs bosons. 

In the second case, the Higgs boson is produced in the 
decay of heavier neutralinos and charginos (xj /X3.4) into 
lighter ones (xt /Xi 2)- To have a large rate from gluino 
and squark decays, the heavy neutralinos and charginos 
need to be gaugino-like whereas the light ones Higgsino- 
like. This corresponds to a /i-term which is small in com- 
parison with the gaugino masses, and to weak gauginos 
which are lighter than squarks and gluinos. 

Fig. [1] shows the typical netralino mass spectrum for 
the "big" and "little" cascade scenarios. 
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FIG. 1: Top: Neutralino mass spectrum in the little cascade 
scenario where Xi,2 are non-degenerate and gaugino like for 
large jj,. Bottom: Neutralino mass spectrum in the big cascade 
scenario where Xi,2 BuTS almost degenerate and Higgsino like 
for small /x. Blue arrows show the decay mode relevant for 
the present search. 

In our analyses, we consider two benchmark SUSY 
models with parameters given in Tab. [l] Benchmark 1 is 
an mSUGRA model with mi/2 = mo = 400 GeV, Aq = 0, 
tan/3 — 10 and sgn(/i) = 1. It corresponds to a model 
where the Higgs is produced in the little cascade X2 



Xih. The dominant production rate is from the pair pro- 
duction of gluinos/squarks, which is given in Tab.m The 
generic final states are (> 2)q+2 [W/Z/h]+E^'^\ As can 
be seen from Tab. |T] Z bosons are much less frequently 
produced compared to Higgs for the reason we have dis- 
cussed previously; on the other hand, W bosons are much 
more frequently produced. In addition, since gluinos are 
flavor blind, they can decay to third-generation squarks. 
This leads to top quarks being produced in the cascade 
decay and further contributes to the large multiplicity of 
W bosons in events in this benchmark model. 

Benchmark 2 features a very small /i-term, correspond- 
ing to a Higgs which is produced in the big cascade. This 
model is the same as SHSP la in [21]. Due to large 
gluino mass, the production is dominated by the process 
pp — >■ qq' . This also indicates that few top quarks will be 
present in these events. In contrast to Benchmark 1, Z 
bosons are almost equally produced in the cascade decay 
as Higgses. 



Model 


1 


2 




940,910 


1000 


mi 


1000 


1000 


m-g 


949 


2036 


1 


163 


138 




306 


-158 


^3 


-518 


306 




535 


625 




305 


148 


m ± 
X2 


534 


625 


tan P 


10 


10 




512 


150 




2.5 pb 


0.41 pb 


BR(qL ^ h) 


30% 


22% 


BR(gL ^ Z) 


3% 


25% 


BR(gi -> W) 


64% 


48% 


a ■ BR(/i) 


0.29 pb 


0.04 pb 


cr ■ BR{h + W/Z) 


0.47 pb 


0.1 pb 


a ■ BR(W/Z) 


1.04 pb 


0.23 pb 



TABLE I: The relevant masses, cross sections and branching 
ratios for the benchmark SUSY models. The spectrum and 
decay branching ratios were calculated using SUSY-HIT [31] . 
a{g,q) are the 2 — >■ 2 LO cross sections involving g and q, 
which were calculated in Pythia. BR(/i), BR(/i -|- W/Z) and 
IiK{W/Z) are the branching ratios for events with at least one 
Higgs boson but no W/Z boson, with both Higgs and W/Z 
bosons, and with at least one W/Z boson but no Higgs boson 
respectively. Masses are given in GeV. 

These two scenarios are, in a sense, orthogonal to each 
other and capture effectively a large portion of the pa- 
rameter space where the relevant branching ratios, bar- 
ring accidental degeneracies in the spectrum, may vary 
only by a factor of 2-3 with respect to the ones appearing 
in Tab. [H 
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IV. SIMULATION AND SUSY CUTS 

Events are generated, showered, and hadronized in 
Pythia 6.4.24 [33] utilizing the "DW tune" [17]. Initial 
and final state radiation as well as multi-particle interac- 
tions are turned on. We included pile-up events assuming 
a luminosity per bunch crossing of 0.05 mb~^. Supersym- 
metric events are isolated by means of cuts which render 
SM backgrounds negligible compared to the SUSY signal. 
We use the cuts in SUSY analyses designed by the CMS 
collaboration for early LHC SUSY searches [Ml US]- First 
we require missing transverse energy E^^^^ > 200 GeV 
and at least three jets with pr > 30 GeV with pseudora- 
pidity I77I < 3. In addition, only events in which the hard- 
est and the second hardest jets have pT > 180, 110 GeV 
respectively are kept. The hardest jet is also required 
to be within the central tracker fiducial volume, i.e. 
I77I < 1.7. Finally, we require Ht > 500 GeV, where 

Ht = J2i=2PT + E^^^^. The jets in the above cuts are 
raw jets clustered using Cambridge/ Aachen (C/A) algo- 
rithm P2] with a cone-size R — 0.5 on calorimeter cells 
with granularity drjx Scf) = 0.1 x 0.1 between —3 < 77 < 3. 
The jet clustering in our analysis is performed using the 
FastJet(v-2.4.2) libraries. 

The preselection cuts and their associated cumulative 
efficiencies on the SUSY signal events from models 1 
and 2 can be seen in Tab. [ll] While relatively robust 
on the SUSY signal events, these cuts are far out on the 
tails of SM QCD, di-boson, and t-tbar backgrounds. In 
the remaining sample, the primary obstruction to recon- 
structing a hadronically decaying Higgs are the SUSY 
events which include W^'s and Z's. To reduce this back- 
ground, we turn to more sophisticated jet substructure 
algorithms. 



cut /sample 


1 


2 


E^'"" > 200GeV 


80.64% 


80.54% 


iV, >3 


75.32% 


78.87% 


Pt,i > 180,PT,2 > 110 


72.29% 


77.72% 


Ht > 500 GeV 


35.54% 


54.47% 



TABLE II: Cumulative efficiencies for the preselection cuts to 
isolate SUSY events. 



V. SUBSTRUCTURE ANALYSIS 

Even in the absence of backgrounds, reconstructing the 
Higgs boson from hadronic jets is generally difficult at 
hadron colliders. However, in supersymmetric produc- 
tion, the Higgs can easily get boosted which collimates all 
of the Higgs decay products into a "fat jet." An effective 
algorithm for identifying these jets and reconstructing 
Higgs candidates is provided by Butterworth, Davison, 
Rubin and Salam (BDRS) p2], where it was used for 
reconstructing SM Higgs. 



This method starts by forming fat jets with cone size 
large enough to capture most of the hadronic products 
of the boosted Higgs; subsequently, one scans within this 
fat jet, looking for a particular jet substructure which 
corresponds to the presumed decay topology of the Higgs 
(in our analysis, — > 2?^ — > Aj). We describe the BDRS 
procedure in detail below: 

(a) Cluster hadronic calorimeter activity into jets by it- 

eratively recombining pairs of closest distance dij. 
For C/A jet algorithm, dij is given by the angular 

distance AR^j = — (f)jY + {Vi ~ ViY- The re- 

combination ends when all objects are separated by 
some minimum ARij > R. Here R is chosen to be 
large enough to contain the decay products of the 
Higgs. 

(b) Uncluster each fat jet into two subjets ji and 

{rriji > mj2). Two criteria must be satisfied by 
these subjets in order to associate the fat jet to 
some presumed heavy parent particle. First, there 
must be a significant mass drop rriji < /i rrij, where 
rrij is the total invariant mass of the parent fat jet 
and /i is a cut parameter. Second, it is required that 
there is no significant asymmetry in the two sub- 
jets defined by: y = min(p|^.^,p|^.J/m2 AR'j^ j^ > 
Vcut- When these two conditions are satisfied we 
exit the loop and dub the jet as candidate Higgs 
jet. 

(c) If the subjets do not satisfy the above requirements, 

then ji is identified as a new fat jet, and step (b) is 
repeated by subdividing ji into a sub-jet pair. This 
is repeated until either prji < 50 GeV or ji can no 
longer be unclustered, at which point the initial fat 
jet is discarded as a candidate for a massive parent. 

In our analyses, we take R — 0.9 — 1.2, /i — 0.5 — 0.667 
and t/cut = (0.3)^, with small variations for different sit- 
uations specified in Tab. |III[ Note, that the BDRS algo- 
rithm is no very effective at selecting the Higgs events: 
only about 10% of the signal survives the initial subject 
analysis. This is due to the fact that the boost of the 
higgs boson is not that incredibly large in these cases. 

Compared to the SM Higgs events, the supersymmet- 
ric events are typically of much greater multiplicity, often 
containing multiple electroweak gauge bosons in addition 
to hard quark and or giuon jets from the decay of squarks 
and gluinos. Many of these hard jets and hadronically de- 
caying W's and Z's can be misidentified as Higgs candi- 
dates, even after the substructure analysis is performed. 
This is the primary obstruction to discovering the Higgs 
as a peak in the jet mass distribution. The issue is com- 
pounded by the fact that the W and Z bosons lie close to 
the mass range expected for the Higgs. Fig. [2] shows the 
typical jet mass distribution obtained for the candidate 
Higgs event. We therefore must supplement the BDRS 
algorithm with additional cuts to suppress the W and Z 
contamination of the Higgs sample as well as other com- 
binatoric jet backgrounds. Different strategies are neces- 
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FIG. 2: The candidate Higgs jet mass distribution for SUSY 
benchmark 1 obtained using the BDRS algorithm. The Higgs 
mass and the t] mass are = 120 GeV and m,, = 10 GeV 
respectively. 



sary for different rj masses. We discuss them in sequence 
below. 

For the low m,, case, the decay products from the rj 
decay are coUimated, and therefore the jet substructure 
is close to the SM Higgs case h — >■ hh. In this case, the 
two subjets from unclustering the fat jet are expected to 
correspond the two ry-jets. In order to reduce the con- 
tamination from W/Z jets, one could consider additional 
cuts on the following variables as discussed in [18]: 

Mass democracy: 

_ min(rnji , mj2) 

omd — r 

max(TOji, TOj2j 

Flow variable: 

Pflow = , , It PT.jS > Pt ■ 

PT,jl+PTj2 

For Higgs decay through two light t^'s, we expect cimd ^ 
1 and /3flow ^ 1- This is based on the fact that both 
higgs and rj are QCD singlets and therefore radiation 
only occurs at the virtuality scale ~ after the rj has 
decayed. The reduced radiation indicates small /3fiow and 
also small shift in the rj jet mass. This is in contrast to 
the QCD jets, where the virtuality scale is governed by 
the initial hard scattering. In [T5] , cuts on these variables 
were used to separate the Higgs jet from the QCD jet. 
In our case, they can instead reduce the combinatoric jet 
backgrounds that are present together with the Higgs. In 
addition, the mass democracy and the flow variable cuts 
are quite useful in distinguishing Higgs jets from W/Z 
jets since the final state radiation in W/Z decay is at a 
larger scale ^ rn^^^z ^ ^ri- For example, in benchmark 
model 1 with m/j = 120 GeV and = 10 GeV they cut 
roughly 75% of the W/Z^s whereas 30% of the Higgses. 

For the high case, the decay products of -q's are less 
coUimated while the two ?7's are more coUimated. This 



makes the four partons more uniformly distributed inside 
the fat jet, giving rise to a truly four-jet decay. This is 
most obvious in the low ruh case, where the allowed phase 
space to decay into 77 is limited. In this case, the two 
subjets found by unclustering the fat jet may not match 
the partonic object from one of the 77's. In addition, due 
to the increased multiplicity of the decay, the subjets are 
typically softer. In order to reduce the W/Z background, 
we need different cuts compared to the light 77 case. 

Number of subjets: The simplest option is to require 
at least four hard subjets inside the fat jet obtained 
from the BDRS procedure: we re-cluster the candi- 
date fat jet into risuhjct subjets with a smaller cone 
size i?sub, 

"■subjet > 4 with Pt > 15 GeV. 

This is easy to understand since W/Z jets typically only 
have two hard subjets. 

Another possibility is to use the planar fiow variable 
introduced in |33j , which is sensitive to whether the color 
flow is linear or isotropic. The planar flow vanishes for 
linear shapes and approaches unity for isotropic depo- 
sitions of energies. In the context of ft. — >■ 277 — >■ 4j, 
the planar flow increases as rriri increases since the flnal 
states become more isotropic. However, in the cases that 
we studied, the number-of-subjet cut is already very ef- 
fective, and we do not include the planar flow in our flnal 
result. 

We have also investigated whether the jet pull vari- 
able [33] significantly enhances signal relative to back- 
ground. We found that, in the cases we analyzed, there 
is little to no improvement as the signal distribution in 
this variable is too similar to the SM gauge boson back- 
ground. However, we have not performed a multivariable 
combined study that could partially enhance the signifi- 
cance [35] . 

In the last step of the reconstruction, a filtration al- 
gorithm cleans up these candidate jets by removing soft 
components. For low 77 mass, one decomposes the fat jet 
to subjets by taking a smaller i?subi and sum up the lead- 
ing Tifiit subjets to obtain the filtered jets. In our analysis, 
we take i?sub = inin(Ai?ji ^2/2, 0.3). For high 77 mass, we 
trim it by only keeping subjets with px > f cut Pt,j [36j . 
where Rgub = 0.2 — 0.3. It should be noted that the 
threshold /cut affects both the accuracy and resolution 
of the Higgs mass. For smaller threshold /cut, mean- 
ing more decay products of the Higgs would be included, 
the reconstructed Higgs mass would be closer to the true 
mass. On the other hand, it is also easier for the con- 
tamination from other softer partons in the same event 
to leak into the Higgs jet, which would worsen the mass 
resolution. The effects of pilc-up events can be seen from 
e.g. Fig. [3| For light rn,^, with pile- up events included, it 
is harderior the fat jet to pass the flow cut. This leads to 
a decrease in the W and Higgs peaks. But on the other 
hand, the continuum background also drops. For heavy 
777,^, there are no qualitative changes in the candidate jet 
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mass distribution. For convenience, we present our final 



result in the figures of Sec. VI without pile-up events 




80 100 120 140 160 
Candidate Higgs Jet Mass (GeV) 



FIG. 3: The candidate Higgs jet mass distribution for SUSY 
benchmark 1 with (solid) and without (dash) pile-up events. 
Top: The Higgs mass and the rj mass are mh = 120 GeV 



and rrirj 
and m„ 



10 GeV respectively. Bottom: — 120 GeV 
30 GeV. The plots are generated using 10 fb~^. 



Substructure cuts are given in Table IHI] 



Tab. mil shows the substructure cuts that we use in 
our search. We will discuss the details in the next sec- 
tion. The concrete values of the cuts vary case by case 
depending on the Higgs and the rj mass, but they are not 
optimized yet. 



mh, nir, 


(120, 10) 


(100,10) 


(120,30) 


(100,30) 


R 


1.2 


1.2 


1.0 


0.9 




0.667 


0.667 


0.667 


0.5 




> 0.7 


> 0.8 


> 0.4 


> 0.4 


/3flow 


< 2% 


< 0.5% 






„min 

Pt 


2.0 


1.0 






Rsub 






0.25 


0.25 


^subjct 






>4 


> 4 


„niin 
PT,sub 






15 


17 



TABLE HI: Jet substructure cuts for different scenarios. Mo- 
menta and masses are in unit of GeV. The cuts in the column 
(100, 10) are only for benchmark 2. 



VI. RESULTS 

We now apply this method of Higgs reconstruction to 
the two SUSY benchmark models for different Higgs and 
rj masses. The Higgs appears as a resonance peak in 
the jet-mass distribution of the fat jets which survive the 
substructure cuts. While the substructure analysis is rea- 
sonably successful at removing hadronically decaying W 
and Z bosons, significant contamination of the sample in 
the 80 — 90 GeV region from these resonances remains. 
The low Higgs mass region, where LEP could have missed 
the Higgs thus remains especially challenging. We con- 
sider separately two different Higgs mass regions: high 
mass {rrih ^ 115 GeV) and low mass {mh < 100 GeV). 

In the heavy Higgs mass region, there is little interfer- 
ence from W and Z contamination of the fat jet sample 
since the peaks in the jet mass distribution are well sepa- 
rated. In this case, one does not need to completely sup- 
press the contribution from hadronically decaying W's 
and Z's, and lower luminosity will be sufficient for Higgs 
discovery. 

In the low mass region, the W and Z jet mass peaks 
share significant overlap with a potential Higgs signal, 
unless the contamination of W and Z bosons can be sig- 
nificantly reduced without losing too much of the Higgs 
signal efficiency. This is in principle possible, due to the 
different decay topology of these events, although issues 
arise when the rj is too light. 

In the case of light Higgs and heavier r], the 4 subjets 
arising from the two rj decays are often resolvable. Ad- 
ditional cuts on the number of sub-jets appearing within 
the fat jet are therefore effective at removing W^'s and 
Z's, even for a relatively light Higgs boson. In the bench- 
mark models we consider, the W and Z background is 
low enough to identify the Higgs. 

For the scenario of both light rj and Higgs mass below 
100 GeV, we find that we cannot remove a large enough 
fraction of the W and Z boson events to be assured that 
an excess in this mass range is due to a Higgs. This is due 
to the fact that light ry's will have substantial relativistic 
boost and correspondingly coUinear decay products. The 
Higgs decay then appears to have di-jet substructure, just 
like the SM gauge bosons. Substructure cuts therefore 
reduce both signal and background to a similar degree. 

One approach to remedy this could be to try to under- 
stand the details of the SUSY background and subtract 
it. This could be done for the Z boson, for example, by 
measuring the number of reconstructed Z bosons in the 
leptonic decay channel. Unfortunately, this can not be 
done with the W boson, since the semi-leptonic decay 
involves a neutrino whose momenta is lost along the LSP 
contribution to the total missing px- Even armed with 
perfect knowledge of these rates, such subtractions are 
especially susceptible to systematic uncertainties in the 
shape of the W and Z jet mass distributions. Due to 
these difficulties, we do not attempt such a subtraction. 

To illustrate the effectiveness of the substructure cuts, 
we looked into our data sample and identified the asso- 
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ciated heavy object for a given candidate fat jet. The 
associated object is defined to be the closest heavy ob- 
ject within an i? = 0.4 cone around the jet. Given that 
information, we are able to count the number of "cor- 
rect" Higgs jets in ±5 GeV window around the true mass, 
and the total number of candidate jets in that window. 
Similarly this can be done for W and Z bosons. These 
numbers can be compared with the number of Higgs, W 
or Z bosons in the sample without subjet cuts to get an 
estimate of the efficiency and the discrimination power. 
As can be seen in Tab. |IV[ a factor of ^ 20 gain in ef- 
ficiency can be achieved for Higgs against W and Z for 
rrijj = 30 GeV, while a factor of ^ 5 for rrir, = 10 GeV. 





Model 1 


Model 2 




(100, 10) 


(100,30) 


(100, 10) 


(100,30) 




before 


after 


before 


after 


before 


after 


before 


after 


H 


6974 


324 
473 


6587 


69 
103 


22450 


700 
831 


22564 


298 
403 


W 


22668 


366 
581 


22435 


7 
26 


63641 


356 
564 


62775 


34 
274 


z 


1296 


18 
390 


1244 




67 


22977 


136 
671 


22933 


19 
269 



TABLE IV: Subjet cut efficiencies for Higgs and W/Z bosons 
in the window ±5 GeV around their true masses. The num- 
ber before cuts are the number of Higgs, W or Z in the event 
sample after applying preselection cuts on the 10^ raw events. 
The number after cuts is presented in the form | where a is 
the number of "correct" Higgs, W or Z jets and b is the to- 
tal number of candidate jets in the respective mass window. 
The "correct" Higgs, W or Z jets are defined as those candi- 
date jets where the closest heavy object within R — 0.4 cone 
around the jet is Higgs, W or Z. 



A. Low T] Mass (m,, = 10 GeV) 

For the low case, we use the modified BDRS 
method with mass democracy cuts and flow cuts to iden- 
tify Higgs jets. As discussed above, the substructure 
analysis is not be able to substantially reduce the con- 
tribution of W and Z bosons while preserving Higgs 
signal events. For both benchmark models with = 
120, 100 GeV, we find candidate Higgs jets and construct 
the jet-mass distribution. 

For benchmark model 1, we take the values for the cut 
parameters to be i? = 1.2, aun > 0.7, /3flow < 2% and 
n-mt — 3. The results for 100, 000 raw events normal- 
ized by the cross section are shown in Figure [4] for both 
high Higgs mass (top panel) and low Higgs mass (bottom 
panel). In this plot, the Higgs mass peaks are well above 
the background and its position is consistent with the 
true Higgs mass. The peaks in the vicinity of 80 GeV are 
from hadronically decaying W's which evade the above 
cuts. To calculate the significance of the Higgs peak, we 
must provide an estimate of the backgrounds from both 
SM and SUSY. The SM backgrounds are negligible as we 
discussed before and are taken to be zero for simplicity. 



while the SUSY backgrounds can be estimated from the 
continuum under the Higgs peak in the jet-mass distri- 
bution. For example, for the case with nih — 120 GeV, 
we take the — 2/-I-1 bins around the peak 120 GeV as 
the signal region and the two adjacent bins for back- 
ground estimation. We find that a Scr discovery of the 
Higgs boson for ^ 10 fb""'^ is possible. For the case of low 
Higgs mass, in the bottom panel of this Figure, these two 
mass peaks are closer. Taking the excess in the ±1 bins 
around the peak 100 GeV as the signal, a 5a significance 
can also be achieved with the same amount of data. For 
an even smaller Higgs mass, the signal peak would be- 
gin to merge with the W peak. Unless the W fake rate 
can be further reduced with additional novel techniques, 
it seems unlikely that a Higgs with mass much smaller 
than 100 GeV can be identified. 



Jet Mass I 
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Candidate Higgs Jet iUlass (GeV) 

FIG. 4: The candidate Higgs jet mass distribution for SUSY 
benchmark 1. Top: nih = 120 GeV and = 10 GeV. 
Bottom: = 100 GeV and = 10 GeV. Events with > 7 
jets {pt > 30 GeV) are vetoed in the bottom plot. 



For benchmark model 2, the results are shown in Fig- 
ure [Sj In the top panel, we use the same cuts as for 
benchmark model 1 and we can see that a 5a discov- 
ery can again be achieved (using — 2/-I-1 bins for sig- 
nal) for ^ 10fb~^ integrated luminosity. In fact, in this 
case the Higgs bosons are generally more boosted due 
to the larger neutralino mass difference. This leads to 
a higher reconstruction efficiency than for benchmark 1, 
and even without the flow cuts we can obtain similar re- 
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suits with smaller luminosity. However, for the low Higgs 
mass, the distribution obtained from using the same cuts 
show a plateau between 80 — 100 GeV. This is due to 
the superposition of W, Z and Higgs contributions. Im- 
posing stronger cuts Q!md > 0.8 and /3flow < 0.5% with 
p™™ = 1 GeV, lead us to the second plot in Figure 
[5j While the W peak is now significantly suppressed, 
and the big peak located around 100 GeV suggests the 
presence of the Higgs boson, the subtraction of the Z- 
background is needed in this case. Naively using the 
same prescription for calculating the significance, we find 
5(7 discovery can be achieved with ^ 25 fb^^ integrated 
luminosity. 




120 140 160 180 200 
Candidate Higgs Jet iVIass (GeV) 



FIG. 5: The candidate Higgs jet mass distribution for SUSY 



benchmark 2. 
Bottom: - 



Top: nih = 120 GeV and 



10 GeV. 



100 GeV and m„ 



10 GeV. 



The resulting candidate Higgs jet-mass distributions 
can be seen in Figures |6|7[ Different from the low rj 
mass cases, the continuum background is small in the 
low mass region and the W/Z peaks are no longer visi- 
ble. This indicates that the cut on the number of sub- 
jets is very efficient in reducing the W/Z contamination. 
But other combinatoric jet configurations can potentially 
leak through the cut since these may have more than two 
hard components and can give rise to a large jet mass. To 
suppress these combinatorics, we use a slightly smaller R 
parameter for the jet clustering algorithm, and include a 
mild cut on the subjet mass democracy omd as shown in 
Tab. 



HI 



For benchmark 1, we require maximum 7 jets 
in the events to further suppress the combinatoric back- 
ground since there are lots of top quarks in the events. 
For the high Higgs mass case, the Higgs peaks are well re- 
constructed, and in both benchmarks a 5a discovery can 
be achieved with roughly 10 and 25fb~^ respectively (us- 
ing — 2/-I-1 and ±2 bins for signals). The results for low 
Higgs mass are similar, but more luminosity (> 35 fb^^) 
is needed due to smaller signal efficiency. 



B. High ri Mass (m^ = 30 GeV) 




120 140 160 180 200 
Candidate Higgs Jet Mass (GeV) 



Moving to the high r/ mass case, the decays of the 
Higgs are more four-jet like. We use the BDRS algorithm 
supplemented with a cut on the number of subjets to find 
the Higgs-like jet. We re-cluster the candidate fat jet into 



subjets using i?sub — 0.25 and require n., 



subjet 



> 4 hard 



subjets with px > 15 GeV. The final candidate Higgs jets 
are obtained after trimming with threshold /cut — 1.5%. 
For the low-mass Higgs, the cuts are slightly adjusted as 
seen in Tab. IhTI 



FIG. 6: Candidate Higgs jet mass distribution for SUSY 



benchmark 1. 
Bottom: mu - 



Top: nih = 120 GeV and 



30 GeV. 



100 GeV and m„ = 30 GeV. Events with > 8 



jets (pt > 30 GeV) are vetoed in both plots. 
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FIG. 7: Candidate Higgs jet mass distribution for SUSY 
benchmark 2. Top: = 120 GeV and = 30 GeV. 

Bottom: = 100 GeV and nir, = 30 GeV. 



VII. CONCLUSIONS 

Purely hadronic Higgs decays of the form h ^ 2r] —>■ Aj 
present a formidable, but not insurmountable challenge 
at the LHC. Such Higgs bosons, when produced in asso- 
ciation with massive exotica, can in many cases be recon- 



structed. In supersymmetric buried Higgs models, cuts 
on events with high pT jets and missing energy remove 
the QCD backgrounds that would otherwise swamp the 
Higgs resonance. 

Once QCD events have been stripped away, the re- 
maining issue is to identify a Higgs resonance among the 
hadronic activity that occurs in typical SUSY cascade 
decay chains. We have applied a substructure analysis 
that aids in the removal of remaining combinatoric back- 
grounds and/or dijet pairs arising from weak gauge boson 
decay. 

For buried Higgs bosons in the range 95 — 120 GeV, 
with both heavy (30 GeV) and light (10 GeV) rj mass, 
discovery at the 5a level is possible with 10 fb^^ of 7 
on 7 TeV LHC running. At lower values of the Higgs 
mass (and low 77 mass), when the Higgs peak overlaps 
significantly with the W and Z resonances, substructure 
cuts do not sufficiently reduce weak gauge boson back- 
grounds. In this case, discovery will be far more chal- 
lenging, and will require more sophisticated techniques 
than those presented here. 
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